Comment on stability of atoms and nuclei in multidimensional TeV scale gravity 
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We discuss stability of atoms and nucleons in the presence of multidimensional gravity 
characterized by the unified energy scale M4+„c 2 R_ 1 TeV. We point out that the multidi- 
mensional gravitational attraction can create additional bound states deep inside of atomic 
O^l ' and nucleon constituents. These states would be characterized by sizes comparable to the TeV 

| scale distance Rt = ft/(M4+ n c). We show that shallow gravity induced bound states between 

(£) . an electron and a quark are ruled out, because collapse processes in atoms to such bound 

states would occur on time scales which are in contradiction to current data on the stability of 
ordinary matter. The gravitational attraction may also produce diquark states, which should 
be considered in addition to those discussed in the literature. The estimates indicate that, for 
TeV scale gravity, the problem of UV divergencies and UV renormalization is crucial. Some 
universal renormalization mechanism should exist, which stabilizes atoms in the presence of 
multidimensional gravity. 
\Q 1 PACS : 04.50.+h 

Recently the idea of multidimensional TeV scale gravity has been put forward in Refs. jj| as a possible 
resolution of the hierarchy problem and as a step toward unification of the Electro- Weak interactions and 
Gravity. One of the main consequences of this approach is that the gravitational potential (fci^) produced by 
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. a point-like mass M will deviate from the Newton's law (Png i r ) — —GM/r very significantly as 

Ph! <Po(r) = -G— ^— j (1) 

\ at distances r <C R n . Here R n stands for the typical size of the n internal dimensions of the (4 + n)— di- 

mensional spacetime and G denotes the observable Newto n's gravitational constant. The value R n is fixed by 
the 4— dimensional Planck mass Mpi = Mu\pi = yKcJG and energy Mpic 2 « 10 19 GeV, and by the energy 
scale i\-4 + „c 2 of the (4 + n)— dimensional gravity. Specifically, for a toroidal internal space 

X 

b: h (Ml , V /n , 17 / 1 TeV \ / Mp, \ 

, R n = -^El- = 2.0 x 10" 17 cm) —P- . (2) 

Eq. (^) indicates that for Mi +n <C Mpi, the size of the compactified dimensions R n becomes much larger 
than the Planck scale Lpi = H/cMpi w 10~ 33 cm. For example, in the case Mi +n c 2 = 1 TeV [Q, one finds 
i?2 ~ 10 _1 cm, i?3 w 10" 6 cm, corresponding to distances below the reach of present experimental tests of 
gravitational forces. 

In this letter, we point out that so large values of R n imply that physics of atoms and nucleons (bound 
to a (3 + 1)— dimensional world brane) becomes very sensitive to the UV behavior of gravity in 4 + n dimen- 
sions. Specifically, atoms and nucleons may become unstable with respect to a collapse. This can be avoided 
by introducing a cut-off for the potential (|l]) below some distance r s , determined by a specific screening or 
renormalization behavior. Below we will restrict the analysis to order of magnitude estimates of the required 
values of r,. 
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We start our consideration by noting that a non-relativistic quantum mechanical system with singular at- 
tractive potential ~ \/r 1+n and n > 1 is known to have no ground state ||. The system is unstable and 
a collapse will occur. The collapse in a relativistic system described by the Dirac equation with potential 
U(r — > 0) ~ l/r 1+ ™ occurs as well jj). Clearly, at atomic and nuclear distances the strength of the additional 
gravitational potential (Q), e.g. between nuclei and shell electrons, is some ten orders weaker than the strength 
of the nuclear Coulomb potential^ ipc = Z\e\/r , e < . Even weaker it is compared with the strength of 
strong interactions between constituent quarks inside the nucleons. Nevertheless, regardless of the strength of 
the l/r 1+n potential (g) at atomic or nuclear distances, a collapse will occur if the form of the potential (|l|) 
holds down to r = without lower bound. 

A possible way to exclude a collapse could consist in some universal renormalization mechanism which makes 
the potential (|lj) non-singular for r < 7' s . For our subsequent considerations we will assume that such a 
renormalization is present. We consider, for simplicity, the following effective potential 

?G(r) = -G^-(^-Y, r s >0, (3) 
r + r s \r + r s ) 

instead of ([!]). 

We assume that the cut-off follows from string theory/M-theory. A natural scale for the upper energy cut-off 
is, then, the string tension M s c 2 . Thus, the distance h/M s c determines the cut-off scale r s . Obviously, such 
a cut-off leads to a lower bound for the discrete part of the energy spectrum. In this regard, several situations 
for the gravity induced bound states are possible: (i) there are "shallow" bound states characterized by energy 
values E = mc 2 — Eb, such that —mc 2 < E < mc 2 , (Eb > denotes the binding energy); (ii) there are deep 
bound states only with E < —mc 2 residing in the positron continuum: (iii) no bound states exist in the potential 

An instability of atoms due to large extra dimensions will occur in the case (i), if the additional gravity 
induced bound states are empty. Then a transition/collapse of an electron, e.g., from the K-shell to such a 
bound state (formed inside the nucleus) can occuiQ. In the case (ii), the bound states will be filled by electrons 
quite rapidly from the vacuum, so that their partners — positrons are emitted to the outside (this process is 
similar to the generation of a charged vacuum in the vicinity of overcritical nuclei f|,|||). As we will see below, 
these bound states are formed at distances where the strength of the gravitational potential overcompensates 
the Coulomb potential. Accordingly, positrons will equally be attracted to practically the same bound states. 
The filling of these positron states from the vacuum will be accompanied by the emission to infinity of their 
partners — electrons. As a result, the system will contain real electron-positron pairs, which are located in the 
vicinity of the gravitating center, so that the total charge of the nucleus does not change. In accordance with 
the Pauli's exclusion principle, no transitions to the filled deep bound states can occur from atomic orbitals, 
so that the stability of atoms is insured in this case. Thus, the stability of atoms is not affected at all by the 
presence of deep bound states. Obviously, no gravity induced transition/collapse occurs in the case (iii). 

In this letter, we consider possible implications of the presence of shallow bound states with —mc 2 < E < mc 2 , 
because, as we have noted above, only in this case the collapse problem becomes relevant. In what follows we 
are confronted with the problem of an adequate description of our system. To be specific, let us consider an 
atom. If the collapse occurs, e.g., for an electron from the K-shell to one of the hypothetical gravity induced 
bound states deep into the nucleus, then the initial state is a well defined atomic state. Due to its smallness 
at atomic distances, the additional gravitational potential can be taken into account as perturbation, and the 



1 An estimate for Hydrogen (Mc 2 = 0.94 GeV, mc 2 = 0.51 MeV ) and the fundamental scale M i+n c 2 = 1 TeV gives, 
e.g., at the Bohr radius r s : ^.^T = (f^)" ~ 1(r? x 4 n 10- 9n and at the proton surface (r p « 10- 13 cm): 

^e>) m « io~ 7 x io- 4n . 

>Pc( r p)e 

2 Due to its small "orbit" size, such a gravity induced bound state should be formed between an electron and a nucleon 
constituent, i.e., a quark. Clearly, the form of the effective gravitating potential will deviate from (^) due to the mass 
distribution inside the nucleus. At atomic distances of order of the Bohr radius the gravitating mass of the center is 
given by the mass of the whole nucleus. Inside one of the nucleons it is given by the mass of one of the constituent 
quarks (e.g., for a u-quark, by m u ~ 310 MeV whereas at smaller distances it will be defined by the bare mass of 
the corresponding quark (m u ~ 4 MeV). This means that a screening form factor should be included into the potential 
(H). For our rough estimates we do not take into account such a form factor. 
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corresponding non-relativistic low-energy physics is well described by the Schrodinger equation. The problem 
arises when one considers the final state — one of the hypothetical gravity induced bound states. We will see 
below, that if such states with energies —mc 2 < E < mc 2 indeed exist, then the corresponding potential energy 
and kinetic energy of the bound electron would be of order of 1 TeV or even higher. Clearly, at such high 
energies the bound state approach is not completely adequate^. We understand that the following estimates 
can only be interpreted as a first crude indication of the collapse problem in TeV scale gravity. 

In traditional models of heterotic string theory, the energy M s c 2 is only few orders less than the Planck energy. 
Accordingly, the cut-off distance r s is only slightly larger than L pi and much smaller than the characteristic TeV 
scale distance Rt- Lp\ ~ 10 _33 cm < r s <C Rt ~ 10 _17 cm. In this case, as will be shown below, the potential 
(||) can support plenty of the deep bound states with E <C —mc 2 . Shallow bound states are possible as well. 
Thus, the cases (i) and (ii) become possible. We should also note that, if the number of deep bound states is 
very large, the gravitational potential may become so strongly modified by the correspondingly large number 
of bound electrons and positrons, that an absolute instability with respect to gaining more and more particles 
from the vacuum develops. The realization of the latter scenario depends significantly on the ratio t s /Rt. In 
this letter, we will not discuss this issue in detail. We only note that, in order to avoid the vacuum instability, 
the scale r s must be much larger than 10 _33 cm. Thus, the string energy must be lowered considerably. Such a 
lowering of the string energy scale M s c 2 down to 1 TeV has been considered, e.g., in HQ. As discussed in |]]], 
only certain types of strings can satisfy this requirement without producing contradictions with observations. 

A natural criterion for the existence of bound states in any attractive potential with a potential energy 
U(r) can be obtained by comparing it with the kinetic energy T^in of the particle confined in a space region 
characterized by some linear dimension r. In general, no relativistic bound state exists if the condition 

|f/(r)|«- (4) 
r 

holds. In TeV scale gravity, this condition (||) is already violated by the unrenormalized potential (|l|) (produced 
by the mass M and attracting the mass m ) at the distance r c given by |Z7(r c )| « hc/r c . Employing 
U (r) = tpG(r)m and excluding R n from ([!]) by means of Eq. (|^) we find 

GmM\ 1/n n / Mm^ 1/n 



The quantity r c represents a rough estimate for the size of the largest "orbit" of a gravity induced bound state 
in the potential (Q). For M 4+n c 2 ~ 1 TeV and Mc 2 = 0.94 GeV, mc 2 = 0.51 MeV, this gives the estimate 
r c w 10 -10 /"i?T ~ 10~ 17 ~ 10 /"cm. From eq. (||) we see that for particles with masses Mc 2 ,mc 2 < M^ n 
the distance r c is smaller than the TeV scale length r c < Rt- This means that the characteristic size of the 
largest " orbit" lies already in the region of strong TeV scale gravity where (multidimensional) quantum gravity 
processes become relevant. These estimates confirm the relativistic ansatz for the kinetic energy T^in ~ hc/r of 
the bound particles which becomes much larger than their rest energies. Accordingly, the rest masses M and m 
in the potential U(r) = to0g(t") should also be replaced by the corresponding expressions for the mass-energies 
fi/(cr) ||, so that effectively U(r) = — (G/r)(h/(cr)) 2 (R n /r) n should be employed for estimating r c . As a 
result, one arrives at the relation r c rj Rt Q for any value of n. 

The above analysis implies that the model is very sensitive to the UV behavior induced by the extra dimen- 
sions. Obviously, for r s 3> r c no bound states will be formed in the rcnormalizcd potential (|3|). In the opposite 
case, when r s «C r c , the potential (Q) should contain plenty of bound states. The situation r s ~ r c is marginal. 
It is, however, obvious that there is a certain critical value r' s — r )r Cl with 7 being some numerical factor of 
order of 1, so that for r s > r' s the bound states disappear, and for r s — * r' s from below the bound states are 
shallow and approach the lower boundary of the upper part of the continuous spectrum. In order to find 7, 
the exact problem should be solved. 



3 More adequate could be a Bethe-Salpeter approach with inclusion of electroweak interactions, and of the Higgs sector 
as well as the of whole tower of massive Kaluza-Klein gravitons. Additionally, at distances r < Rt and energies above 
M4+„c 2 , strong TeV scale quantum gravity effects, which could lead to topology fluctuations, should also be taken into 
account. 

4 In this case, the condition |C/(r c )| « hc/r c can be written as ~4>g{t c ) /c 2 w 1 which implies that the gravitational 
potential is in the range of intermediate gravity and at the upper bound of the employed Newtonian approximation. 
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We note that the existence of empty shallow bound states would strongly destabilize ordinary atoms on a 
time scale which is comparable with times of usual atomic transitions. Indeed, let us assume that there is a 
TeV scale gravity induced shallow bound state ipG( r ) for an electron. In rough approximation, such a state 
can be modeled by the Dirac equation with an effective potential corresponding to a spherical well of a depth 
Uq rs ra|<^>G(0)| and of the radius r s . This potential should transform into the standard Coulomb potential at 
r > r ~ r c Q It is important that, regardless of the details of the problem on distances r ~ r s , the behavior 
of ipG( r ) outside the gravitational potential is determined by the Coulomb potential and by the value of the 
energy E. For shallow bound states the long distance asymptotic is given by 

— \r 

iP G (r) = A - , (6) 



where A = ^/m 2 c 2 — E 2 /c 2 /h jj); Ao stands for some constant which should be determined from the solution 
at small distances. To be more specific, the solution in the region where gravity dominates ( r ~ Rt ) must 
be matched with the standard behavior in the Coulomb potential || . In some sense, the constant Aq defines 
the part of the gravity induced bound state which spreads outside the region of strong gravitational attraction. 
The quantity A -1 sets a typical distance scale for the extension of V'G( r )- Obviously, A -1 3> r c rs Rt for 
an electron. The amplitude Aq in the asymptote (0) determines also the scattering phases ||, and for any 
realistic potential these phases should be finite (see M), and therefore Aq should also be finite in the limit 
r s -► 0. 

The above analysis implies that the overlap of the gravity induced shallow bound state with any atomic state 
(with characteristic extension r in rs 10 _8 cm) is determined by the ratio \ = ( r m^) -3 > 10~ 8 rather than by 
the smallness (r c /ri n ) 3 ~ 10 -27 . For the K-shells of heavy elements, the ratio \ can be as large as 10~ 2 . 
Thus, the corresponding matrix elements of the photon assisted transitions from atomic orbitals to gravity 
induced shallow bound states are only insignificantly reduced by x- This becomes especially clear when one 
considers observational implications of such transitions/collapses. Indeed, in order to reconcile the existence of 
the shallow bound states with current data on the stability of ordinary matter, the collapse rate must be smaller 
than at least the inverse age of the Universe 10 _18 s _1 . In fact, observational limitations on the stability are 
much stronger. Experimental lower bounds on a possible disappearance of K-shell electrons in Iodine atoms 
yield, e.g., a half-life time of such processes > 10 29 — 10 31 s [[}). The above value of \ can reduce the typical 
rates of the atomic transitions from 10 15 — 10 8 s _1 to about 1 — 10 7 s _1 only. This, implies that the existence 
of shallow gravity induced bound states (the possibility (i)) is absolutely excluded for electrons. 

Thus, if multidimensional TeV scale gravity induces bound states, then these states should not be shallow. 
In fact, these states should lie deep enough so that the escaping part of the electron wave function (^) does 
not significantly modify the probability of finding the atomic electrons at the nucleus. Otherwise, such deep 
gravity induced states would contradict to the current data on the atomic hyperfine splitting, which is known 
to be sensitive to the probability of finding electron near the nucleus H. In this letter, we will not further 
discuss possible implications of the hypothetical gravity induced deep bound states. We only note that in the 
case (ii), the Pauli exclusion principle would stabilize the system, so that transitions of atomic electrons to such 
states will not occur. Thus, in the case (ii), the presence of the TeV scale gravity potential will have no direct 
consequence on the stability of atoms. 

We note that, for a realization of the case (iii) with no bound states present at all, the cut-off r s should 
satisfy the condition 

r s > 7 r c rs R T . (7) 

The above analysis has been carried out for electrons in the field of some massive center. A similar stability 
problem, which, however, leads to less drastic consequences than in the case of the atomic electrons, exists for 
quarks in nucleons. In this case, the TeV scale gravity may induce, e.g., diquark states at various energy scales. 
Such states, which are characterized by the quarks separation ~ Rt should be considered in addition to those 
discussed in the literature pOl. We will not dwell on this here. 



5 The distance ro, where the potential energy of an electron in the gravitational field of another particle becomes 
comparable with its energy due to the Coulomb attraction by a unit electric charge, is in the ultra-relativistic limit 
roughly given by (G/ro)(h/(cro)) 2 (R n /ro) n ~ e 2 /ro (without running couplings taken into account). Accordingly, one 
has r w {hc/e 2 ) 1/{ - n+2) R T = (137) 1/(n+2) 7? T > r c rs R t . 
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Summarizing the above discussion, we are led to the conclusion that the UV problem is very acute in multi- 
dimensional TeV scale gravity^, and any renormalization scheme should ensure that at least no shallow bound 
states with binding energies below the electron rest mass can be formed. A similar stability problem as for TeV 
scale gravity will also occur in Randall- Sundrum models p[ . Although the effective on-brane gravitational 
potential on the Planck brane differs from that on a probe brane, the short distance behavior on both types 
of branes is given by a r 1+n term with n > 1 |^2| , so that an appropriate circumvention of the corresponding 
collapse problem should be found. We should also note that our results are not valid for models |0|], where 
due to the induced metric on the brane, the gravity on the brane is exactly 4— dimensional all the way down to 
distances of order h/Mpic. 

Acknowledgments We thank Valery Rubakov and Gia Dvali for valuable comments. This work has been 
initiated during A.Z.'s visit at the College of Staten Island of CUNY. A.Z. acknowledges support by CUNY, 
the US Federal Department of Foreign Affairs, the programme SCOPES (Scientific co-operation between East- 
ern Europe and Switzerland) of the Swiss National Science Foundation, project No. 7SUPJ062239, and the 
Spanish Ministry of Education, Culture and Sport. U.G. acknowledges financial support from DFG grant KON 
1575/1999/GU522/1. 



272, 



[2] 
[3] 

[4] 
[5] 
[6] 
[7] 



[8] 
[9] 
[10] 

[11] 
[12] 
[13] 



hep-ph/9803315 



N. Arkani-Hamed, S. Dimopoulos and G. Dvali, Phys. Lett. B429, (1998), 263 
N. Arkani-Hamed, S. Dimopoulos and G. Dvali, Phys. Rev. D59:086004 (1999), |hep-ph/9807344 ~ 
LD. Landau and E.M. Lifshitz, Quantum Mechanics, v. Ill, Pergamon Press, Oxford-NY, 1977. 
V.B. Berestetskii, E.M. Lifshitz and L.P. Pitaevskii, Relativistic Quantum Theory, Part I, v. IV, Pergamon Press, 
Oxford-NY, 1971. 

M. Lavelle and D. McMullan, Phys. Rep. 279 (1997) 1 
Ya.B. Zeldovich and V.S. Popov, Uspehi Fizich. Nauk, 105 (1971) 403. 

W. Greiner, B. Miiller and J. Rafelski, Quantum e lectrodynamics o f strong fields, Springer, Berlin, 1985. 
JD. Lykken, Phys. Rev. D54 (1996) 3693 - 3697, jhep-th/9603l"33| ; 
I. Antoniadis and K. Benakli, Int. J. Mod. Phys. A 15 (2000) 4237 - 4286, 



65, 



hep-ph/9509344. 



hep-ph/0007226 



C. Misner, K. Thorne and J. Wheeler, Gravitation, Chapter 23, Freeman, San Francisco, 1973. 
Y. Aharonov et al, Phys. Rev. D52 (1995) 3785 - 3792; P. Belli et al, Phys. Lett. B460 (1999) 236 - 241. 
M. Anselmino, E. Predazzi, S. Ekelin, S. Frederiksson, and D. Lichten berg, Rev. Mod. P hys. 65 (1993) 1199 
L. Randall an d R. Sundrum, Phys. 
4690 - 4693, ^ep-th/9906064 . 



1233. 

Rev. Lett. 83 (1999) 3370 - 3373, |hep-ph/9905224 Phys. Rev. Lett. 83 (1999) 



hep-th/0002091; D. Chung, L. Ev erett and H 



Giddings, E. Katz and L. Randall, JHEP 0003 (2000) 023 
Davoudias, "Experimental probes of the Randall-Sundrum infinite extra dimensio n", |hep-ph/0010103 
G. Dvali and G. Gabadadze, "Gravity on a brane in infinite-volume extra space", hep-th/0008054. 



Similar stability arguments for atoms can be applied to the physics on a brane of some finite width Rb- The traditional 
attractive Coulomb potential should convert into its multidimensional counterpart Vc(r) ~ l/r n+1 at r < Rb- 
Accordingly, additional induced states with r c ~ Rb could exist unless some universal renormalization mechanism in 
4 + n dimensions makes the potential less singular at such distances. 
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